Gamma-ray induced cascades and magnetic fields in intergalactic medium 
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We present the results of Monte-Carlo simulations of three-dimensional electromagnetic cascade 
initiated by interactions of the multi-TeV 7-rays with the cosmological infrared/optical photon 
background in the intergalactic medium. Secondary electrons in the cascade are deflected by the 
intergalactic magnetic fields before they scatter on CMB photons. This leads to extended 0.1-10 
degree scale emission at multi-GeV and TeV energies around extragalactic sources of very-high- 
energy 7-rays. The morphology of the extended emission depends, in general, on the properties of 
magnetic fields in the intergalactic medium. Using Monte-Carlo simulated data sets, we demonstrate 
that the decrease of the size of extended source with the increase of energy allows to measure weak 
magnetic fields with magnitudes in the range from < 10^^^ G to lO^^'^C if they exist in the voids 
of the Large Scale Structure. 
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I. INTRODUCTION 

Spectral and timing properties of astronomical sources 
of very high-energy 7-rays could be strongly affected by 
the development of electromagnetic cascades on the way 
from the source to the Earth. These cascades could be 
initiated by interactions of the 7-rays with the ambient 
radiation fields inside the 7-ray source, in the source host 
galaxy and in the Milky Way galaxy, as well as with cos- 
mological photon fields in the intergalactic space. The 
ubiquity of particle cascades has two-fold consequences. 
On one, pessimistic, side, they complicate the interpre- 
tation of the observational data in the very-high-energy 
(VHE, 0.1-10 TeV) 7-ray band. On the other, optimistic, 
side, with enough spectral and angular resolution, one 
can quantify the influence of the cascade on the observed 
source signal and not only reconstruct the spectrum of 
the primary 7-ray source, but also use the information 
about properties of the cascade to study the physical 
characteristics of the medium in which the cascade has 
developed. 

An illustration of such possibility is given by the ob- 
servations of the effect of absorption of VHE 7-rays on 
the extragalactic background light (EBL). Observations 
of this effect in the spectra of distant blazars are now rou- 
tinely used to constrain the EBL spectrum and, in this 
way, the models of cosmological evolution of galaxies and 
stellar populations (see [l| for a recent review). 

In the derivation of the constraints on the EBL from 
the observations of absorption of VHE 7-ray flux from 
distant blazars one usually adopts an assumption that 



secondary emission from the e^e^ pairs, deposited in 
the intergalactic space by the absorbed 7-rays, is not 
detectable. This assumption is valid if the trajectories 
of the secondary pairs are significantly deflected by the 
magnetic fields during their radiative cooling. This is 
true if the strength of the extragalactic magnetic fields 
(EGMF) is higher than - lO^^^ Q [3,[|. 

However, the strength of EGMF is not measured (see 
3 for a review of measurements of cosmic magnetic 
fields), so that there is no direct way to check the validity 
of the the assumption of the non-detection of the cascade 
emission. It is also not possible to constrain the values 
of EGMF from the cosmological models of the origin of 
magnetic fields in the Universe. 

The problem of the origin of 1 — 10 /.jG magnetic fields 
in galaxies and galaxy clusters is one of the long stand- 
ing problems of astrophysics/cosmology (see 0, [a, for 
reviews). On the galaxy scales, such fields are thought to 
be produced via dynamo mechanism [^, from "seed" 
primordial fields of unknown origin. Recent observations 
of the presence of strong, ~ 10 fiG, magnetic fields in 
the high-redshift galaxies [13, [Hi point to the fact that 
the dynamo mechanism is extremely efficient and/or the 
seed magnetic fields are quite strong so that the 10 ^G 
magnetic fields are generated already at large redshifts. 

The dynamo mechanism could work also on the scales 
of galaxy clusters [l3| , although the seed magnetic fields 
would be amplified by a much smaller factor in this case. 
The fact that the observed strength of magnetic field in 
galaxy clusters is comparable to that of the individual 
galaxies [l^ implies that the seed fields in the galaxy 
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clusters should be orders of magnitude higher than the 
seed magnetic fields for the galactic dynamos. Recent 
simulations show that such stronger seed magnetic fields 
could be produced via magnetized winds from the cluster 
galaxies [11, [Hi. 

The origin of the seed magnetic fields for the galactic 
dynamos is largely uncertain. All the models of produc- 
tion of such seed fields invoke the difference of mobility 
of electrons and protons, which results in production of 
electric currents and of the " battery- like" effects at dif- 
ferent stages of cosmological evolution, from the phase 
transitions in the early Universe, to the formation of first 
galaxies [3| ■ A common feature of all the theoretical 
models is a prediction that a larger or smaller fraction of 
the space outside galaxies, galaxy clusters and filaments 
of the large scale structure has to be filled with very weak 
magnetic field, with the strength approaching that of the 
primordial seed fields. Depending on the model, the pre- 
dictions of the typical strength of the EGMF in the voids 
of the large scale structure ran ge f rom < lO^^" G up to 
- 10-9 G (see e.g. [1, [3, mm 111)- Moderate ob- 
servational limit on the present day strength of EGMF 
B < IQ-^G comes from the limit on rotation measure 
of emission from distant quasars [^ . Similar restriction 
comes from the analysis of anisotropics of the cosmic mi- 
crowave background [H, H^] . 

Gamma-ray observations could, in principle, be used 
to constrain properties of the EGMF using the imaging 
[^ and/or timing [U, [2^ of the 7-ray signal. The idea 
is to observe or, at least, constrain the properties of the 
inverse Compton (IC) emission produced by the e'^e" 
pairs deposited by the absorbed 7-rays in the intergalac- 
tic space. If the EGMF strength is below 10-^^ G (plausi- 
ble assumption, in the view of the cosmological models of 
the origin of magnetic fields) , one expects that the defiec- 
tions of the secondary electrons and positrons by EGMF 
should result in the appearance of an EGMF-dependent 
extended emission around initially point sources and/or 
in an EGMF-depcndent time delay of the emission from 
the secondary pairs. If the EGMF is stronger than 
10-1^ G, VHE 7-ray observations still can be used to de- 
rive a lower limit on the EGMF strength, because strong 
EGMF > 10-1^ G can stiU be revealed via the extended 
emission around VHE sources [2| with properties not di- 
rectly dependent on the EGMF structure. 

In what follows we explore the possibility to measure 
the EGMF using the data of VHE 7-ray observations of 
extragalactic sources. In the Section |IT] we extend ideas 
of ref. [^ on dcvelopement of electromagnetic cascade in 
the intergalactic medium and find the range of EGMF 
strengths which lead to observable extended emission 
around extragalactic 7-ray point sources at GeV and TeV 
energies. In the Section Hill we model the development of 
the 7-ray induced cascades in the intergalactic medium 
via Monte-Carlo simulations. This allows us to simulate 
the cascade-induced extended emission around initially 
point 7-ray sources and to study the dependence of the 
properties of this extended emission on the strength of 



EGMF. Using the simulated data sets, we demonstrate 
that the assumed strength of EGMF can be calculated 
from the measurement of the energy dependence of the 
size of the extended emission. We summarize our results 
in the Section HVl 



II. 



ELECTROMAGNETIC CASCADE IN 
INTERGALACTIC MEDIUM 



Multi-TcV 7-rays emitted by distant point sources are 
not able to propagate over large distances because of the 
absorption in interactions with EBL. The pair production 
on EBL reduces the fiux of 7-rays from the source by 



(1) 



where F{E^^) is the detected spectrum, Fq{E^^) is the 
initial spectrum of the source and T{Ejg, z) is the optical 
depth with respect to the pair production on EBL, which 
is a function of the primary photon energy E^^ and of the 
redshift of the source z. 

The e^e~ pairs of the energy E^ produced in interac- 
tions of multi-TeV 7-rays with EBL photons produce sec- 
ondary 7-rays via IC scattering of the Cosmic Microwave 
Background (CMB) photons to the energies 
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where ecMB = 6 x 10 eV is the typical energy of CMB 
photons. In the above equation we have assumed that 



the energy of primary 7-ray is Ej 



2Ee- Upscat- 



tcring of the infrared/optical background photons gives 
sub-dominant contribution to the IC scattering spectrum 
because the energy density of CMB is much higher than 
the density of the infrared/optical background. 

Deflections of e^g- pairs produced by the 7-rays, 
which were initially emitted slightly away from the ob- 
server, could lead to "redirection" of the secondary cas- 
cade photons toward the observer. This effect leads to 
the appearance of extended emission around an initially 
point source of 7-rays. 

The cascade electrons loose their energy via IC scat- 
tering of the CMB photons within the distance 
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Comparing this to the Larmor radius in the magnetic 
field S, 
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one can find typical deflection angle of the cascade pho- 
tons 
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A simple geometrical calculation Q shows that the the 
bulk of the IC upscattered cascade 7-rays of the energy 
arrives within an angle 
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where D is the distance to the source and D^^ is the 
mean free path of the primary 7-rays. The above expres- 
sion is valid for t{E^q,z) > 1. If the extension of the 
cascade- induced "glow" around initially point source is 
determined by the deflections of the cascade electrons by 
the EGMF, the extension of the "glow", Qcxt is expected 
to shrink inversely proportional to the 7-ray energy. 



If t{E^ 



> 1, most of the power of the primary 



photon beam at the energy E^g is transfered to the power 
of electromagnetic cascade, so that the the total flux of 
the extended source is of the order of the primary source 
flux at the energy E^„. The extended source flux at a 
given energy E^ is further modified by the absorption on 
the way toward the Earth, so that the resulting flux of 
the extended source at the energy E^ is Fcxt{Ef, Qext) — 

The above estimate of the flux of the extended source is 
obtained under assumption that the magnetic field of the 
strength B fills most of the intcrgalactic medium between 
the source and observer on the Earth. On the other hand, 
the line of sight can be intersected by the elements of 
the large scale structure, in which the magnetic field is 
much stronger than the typical EGMF. Strong defiections 
of electrons in the parts of the cascade traversing the 
large scale structure could lead to the suppression of the 
secondary cascade flux in the direction of observer. To 
take into account this effect, it is convenient to introduce 
the "volume filling factor" , V{B), of the field of particular 
strength B. The fraction of the path along the line of 

1 /3 

sight, occupied by the field of the strength B is [V{B)] ' 
and the 
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Combining Eqs. ([S]), ([7]), one can sec that a measure- 
ment of the total flux and of the energy-dependent size 
of the cascade emission around an extragalactic point 
source enables to extract the information about the pa- 
rameters of extragalactic magnetic fleld, B, V{B), along 
the line of sight. 

It is easier to detect an extended emission around a 
given point source if the angular size of the extended 
source is not too large. The energy of the detected 
7-rays, E^, and the EGMF, B, enter in the Eq. ^ 
for the source size through a combination E^/B. This 
means that stronger magnetic flelds could be measured 
by the instruments sensitive at higher energies. At the 
highest energies, the sensitivity of a telescope is deter- 
mined mostly by the effective collection area Aas{E^). 
Let us consider, as a reference, a source at the redshift 
z ~ 0.03, producing a flux F{E^g) ^ IQ-^^ erg/cm^s. If 



E^g ^ 100 TeV, the optical depth T{Ejg,z) ~ 10, so that 
all the primary 7-ray power is transfered to the extended 
source emission at the energy ~ 6 TeV. The statistics 
of the extended emission signal is 
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The statistics sufficient for the extended source analy- 
sis is thus achieved within the exposure time Tcxp ~ 
50 hr (typical for the current generation ground-based 
7-ray telescopes) with a telescope with an effective area 
Acff ^ 10 km^. Such effective area is achievable with 
the planned next-generation Cherenkov telescope array, 
optimized for the observations in the 10 TeV energy 
band [1^. Otherwise, a comparable signal statistics is 
achievable with a much smaller collection area, A^s ~ 
0.1 km^, provided that the exposure time is Texp ~ 1 yr. 
This is achievable with the next-generation HAWC water 
Cherenkov telescope HI] . 

Extremely weak magnetic flelds B < 10^^^ G could be 
detectable through the degree-scale extended emission in 
the GeV, rather than TeV energy band, see Eq.®. The 
primary photons which lead to production of extended 
emission in the GeV band have energies E^g ~ 1 — 3 TeV, 
see Eq. ([2]) . The mean free path of such photons is much 
larger than that of the ~ 100 TeV photons, so that de- 
tectable extended emission is expected around distant 
(z 3> 0.1), rather than nearby, sources in this case. Even 
smaller fields can be probed with timing analysis of ob- 
served signal, as was discussed in the ref. [22]. 

The extended emission in the GeV energy band has 
to be detected on top of the diffuse 7-ray background 
(DGRB). The characteristics of this background are un- 
certain a.t p resent. Previous measurements done by 
EGRET [11 ill indicate that a significant part of DGRB 
can be due to unresolved point sources (e.g. blazars). 
The true DGRB should contain at least contribution from 
cosmic ray interactions with CMB, but this contribution 
can be significantly below the EGRET value [11]. In 
near future Fermi [27| will provide a new measurement of 
DGRB. Assuming that the intensity of the DGRB mea- 
sured by Fermi will be by a factor f Fermi < 1 lower than 
the EGRET measurement, we can estimate the sensitiv- 
ity of Fermi for the detection of extended emission around 
high-redshift sources in the GeV energy band. The num- 
ber of DGRB 7-rays within a 9 1 degree region on the 
sky is 
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Re-scaling the Eq. ([H]) for a brighter blazar F{E^g) ~ 
10^11 crg/cm^s one finds that the signal statistics is in 
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this case: 



III. NUMERICAL SIMULATIONS 
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Assuming the exposure time Texp ~ 1 yr, one finds that 
a moderate collection area Aeff < 1 m^, typical for the 
space-born 7-ray telescopes, like Fermi, is sufficient for 
the study of expected extended emission. 

Assuming that the sensitivity of a telescope is sufficient 
for the detection of extended emission, one can estimate 
the maximal measurable magnetic field from the condi- 
tion that the source extension Qcxt{Ej) does not exceed 
the size of the telescope's Field of View (FoV) 0foV- 
Taking the source size at the energy ~ 1 TeV as a 
reference, and substituting OpoV at the place of Oext in 
Eq. ([6]), one can find that the maximal measurable field 
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An increase of the size of the FoV, expected with next- 
generation ground based 7-ray telescopes, such as the 
Cherenkov telescope array optimized for the 10 TcV 
energy band [23j . or the next-generation wide field 
Cherenkov telescope array AGIS [20|, will provide, apart 
from a better sky coverage, an improvement of sensitivity 
for the detection of stronger EGMF. 

The weakest magnetic fields which can be probed is 
found from observation that 0oxt can not be measured 
if it becomes smaller than the size of the point spread 
function of the telescope, OpsF- Since 6oxt ^ (see 
([6])), the largest source extension is achieved at lowest en- 
ergies. Assuming Spsp ~ 0.1° (typical for the Cherenkov 
telescope arrays and for Fermi), one finds from ^ 
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The characteristics of a set-up optimized for the detec- 
tion of weakest EGMF B ~ lO^^*^ G, should be quite 
different from the one optimized for the detection of 
stronger EGMF, B - 10"^^ G. In this case, high sen- 
sitivity at low energies E^ < 0.1 TeV and good angu- 
lar resolution arc required. These can be achieved with 
the planned next-generation low-energy extensions of the 
ground-based Cherenkov telescope arrays, like HESS-II 
[33, MAGIC-II [U and CTA ^ or with Fermi 
Another important difference between the approaches for 
detection of strong and weak EGMF is different choice of 
sources. Strong EGMF could be detected via a study of 
extended emission around nearby {z < 0.1) extragalactic 
TeV sources, while weak EGMF could be detected via 
studies of extended emission around sources of sub- TcV 
7-rays in the distant Universe (z > 0.1). 



The EGMF-determined energy dependence of the size 
of the extended emission in the TeV energy band, given 
by Eq. ([6]), provides a possibility to extract the infor- 
mation on the properties of the EGMF from the obser- 
vational data. In the qualitative discussion of the previ- 
ous section we have assumed, for the sake of simplicity, 
that each primary 7-ray produces electron and positron 
of energies equal to one-half of the primary 7-ray en- 
ergy, Ee = Ery/2, and that electron and positron emit 
IC photons of energies equal to (£'e/?Tie)^ecMB- In other 
words, we have adopted a "monochromatic approxima- 
tion" for the production spectra of e+e^ pairs and of the 
IC scattered photons. A scatter of the energies of cascade 
particles can, in principle, lead to "blurring" of the en- 
ergy dependence of the surface brightness profiles of the 
extended source, thus preventing the possibility of mea- 
surement of magnetic field via the study of the source ex- 
tension. Besides, both spectrum and morphology of the 
extended emission could depend on the spectrum of the 
primary source. This dependence could lead to the loss 
of the characteristic E~^ behaviour of the source exten- 
sion and to the loss of information about the properties 
of EGMF. 

In order to study these problems in more details, we 
have developed a Monte-Carlo code for the modeling of 7- 
ray-induced electromagnetic cascades in the intergalactic 
space, with account of the EGMF. 



A. Primary 7-rays 

We consider, as a reference case, an extragalactic 
source that injects 7-rays with a powcrlaw energy dis- 
tribution, dN^g/dE^g ^ E~^ cxp{—E^g/Ecut) with the 
powerlaw index F = 2 and cut-off energy Ecut = 300 TeV. 
The calculations presented below can be generalized (re- 
peated) for an arbitrary assumptions about the primary 
source spectrum. As it is explained below, an analysis of 
the real data sets aimed at the measurement of EGMF 
has to include Monte-Carlo simulations of the properties 
of the extended sources produced by the primary point 
sources with different spectra. 

To calculate the pair production by the primary 7-rays 
in the intergalactic space, we assume the energy depen- 
dence of the mean free path of the 7-rays calculated in 
the Ref. [s^ for "nominal" Cosmic-Infrared-Background 
(CIB) model. We take into account the uncertainty of 
the measurements of CIB by allowing for different over- 
all normalizations of the CIB. The pair production on the 
CMB becomes important at the energies above ~ 50 TeV. 
Since our calculations extend to the primary 7-rays with 
energies above 50 TeV we also take into account the CMB 
photon background. 

If the mean free path of the 7-rays of the energy Ejg 
is Djg, the probability density of creation e~^e~ pair at a 
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distance x from source is 



C. Extragalactic Magnetic Field. 
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so that the distances Xi at which i-th 7-ray creates an 
e^e~ pair could be expressed through a random number 
Pi distributed in the interval < < 1 as 

x,^ ~D^,-\n{P,). (14) 

Using the production spectrum of the e"*"e~ pairs, 
dN{Ee)/dEf, — f{ujQ,Eyg,Ee), where luq is the energy 
of the background photon, we calculate the energy of 
electron, produced by the z-th primary 7-ray, taking into 
account that the probability for an electron to have en- 
ergy Ee,i is 



(15) 



Having calculated the electron energy E^^i we calculate 
the positron energy as E^+^i = E^^ — E^^i. Initial di- 
rections of motion of electron and positron Vg± coincide 
with the primary photon direction. 



B. 



pairs 



We model the propagation of both electrons and 
positrons in the EGMF by solving the equations of mo- 
tion. 



d"Ve 

leme—— = ±e Ve X B 
dt 



(16) 



here 7e is Lorentz factor of e~ or e"^, me,e are, respec- 
tively, electron mass and charge and Vg is its velocity. 

The IC scattering of the background photons is treated 
in a way similar to the one used for the pair production, 
via calculation of the distance of emission of each next 
IC 7-ray, using equation similar to ()14p with a substi- 
tution Die = (ct'^cmb)""'^, where gt is the Thomson 
cross-section, at the place of D^^. The energy of each 
subsequent upscattered photon is calculated based on 
the differential cross-section of IC scattering 3J], in a 



way similar to the calculation of the energies of the e^e" 
pairs, described above. Having calculated the energy of 
an IC 7-ray, E'^, we decrease the electron/positron en- 
ergy, Ee Ee ~ Ej. We repeat the generation of the 
IC 7-rays for all branches of the cascade until energies of 
electrons/positrons become less then 0.1 TeV. 

The secondary cascade 7-rays produced via IC scat- 
tering process by the e+e~ paris could themselves be ab- 
sorbed by interactions with CIB photons and inject addi- 
tional e"*'e~ pairs. We take this possibility into account 
in our calculations. We do not include the cosmological 
evolution of the CMB and CIB, because we consider only 
the relatively nearby extragalactic 7-ray sources, situated 
at the distances ~ 100 Mpc. 



The energy attenuation length of the e+e^ pairs in the 
intergalactic medium is short compared to the coherence 
length of the EGMF, Degmf ~ 1 Mpc, estimated to be 
of the order of the typical distance between the galax- 
ies [1, [s^ . This means that one can assume that elec- 
trons and positrons propagate in an ordered magnetic 
field (rather than diffuse in a random field) during their 
radiative cooling. In our modeling we divide all the space 
into grid of cubic cells with the side length 1 Mpc. In each 
cell we choose a constant magnetic field with random ori- 
entation. We assume, for simplicity, that the strength of 
EGMF is the same in each cell and leave the investigation 
of the effect of variations of the EGMF strength along the 
lines of sight toward individual extragalactic sources for 
the future work. 



D. Properties of 7-ray-induced cascades in 
intergalactic space 

The results of typical runs of our code are presented 
in the Figs. [T]and[21 Fig. [T] compares the cascades in- 
duced by the 7-rays of the energy E^g — 300 TeV emitted 
from a source at the distance D = 120 Mpc. The three 
cases shown with black color in the figure correspond to 
different choices of the strength of EGMF. For compar- 
ison, we show in the same figure in red a cone with an 
opening angle 6jet = 5° which is believed to be a typical 
opening angle of the blazar jets (6jot — l/rjot, where 
Fjot ~ 10 is the bulk Lorentz factor of the jet). The 
leftmost black-colored cascade develops in the EGMF of 
the strength B = 10^^'^ G, the cascade in the middle 
develops in EGMF B ~ 10~^'^ G while the cascade on 
the right develops in a magnetic field B = 10~^^ G. One 
can clearly see that the increase of the strength of the 
magnetic field leads to the increase of the opening angle 
of the cascade. Also the opening angle of the cascade 
becomes larger than the typical jet opening angle at a 
certain distance from the source. This distance depends 
on the EGMF strength. 

Fig. [2] shows the dependence of the geometrical prop- 
erties of the cascade on the energy of the primary 7-ray. 
The left, middle, and right cascades are initiated by the 
7-rays of the energies 10, 30 and 100 TeV, respectively. 
The decrease of the mean free path of the primary 7-ray 
with the increase of energy is evident from this figure. 
The decrease of the mean free path leads to the earlier 
on-set of the electromagnetic cascade. 



E. Extended emission around an extragalactic 
point source 

To model the detection of the 7-ray-induced cascades 
by the space- or ground-based 7-ray telescope, we adopt 
the following procedure, illustrated by Fig. [31 
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FIG. 1: Examples of showers from 7-rays with primary en- 
ergy Eo = 300 TeV developing in the EGMF of the strength 
10"" (left), 10"^^ (middle) and 10"^^ G (right). Only pho- 
ton tracks with energies above 0.1 TeV are shown. The mean 
free paths of the electrons and positrons are much shorter 
(~ 1 kpc) than the cascade development scale and they are 
not visible in the Figure. 




FIG. 2: Examples of showers from 7-photons with different 
primary energies 10 (left), 30 (middle) and 100 TeV (right) 
in EGMF 10"^^ G. 



Simulations of electromagnetic cascades initiated by 
the primary 7-rays injected by an isotropically emitting 
source or by a source emitting into a jet with a certain 
opening angle take a lot of computing time. Besides, 
most of the cascade photons initially emitted at large an- 
gles with respect to the line of sight would never reach the 
detector, so that the machine time spent on their calcu- 
lation would be wasted. Taking this fact into account, we 
inject all the primary photons in the same direction. We 
save the coordinates and arrival directions of all the cas- 
cade photons with energies above 0.1 TeV, which crossed 




FIG. 3: Scheme of Monte-Carlo calculations of the extended 
emission. 

the sphere with centre at the location of the source and 
with radius D equal to the distance to the detector. 

Next, to take into account the characteristics of the 
detector, we choose a circle of the angular radius dpsF = 
DOpsp where 6'psF is the angular resolution of the detec- 
tor. To calculate the extended emission produced by an 
isotropically emitting source, we place the detector cir- 
cle at a set of random positions on the sphere and sum 
all the cascade photons which pass through the detec- 
tor circle at a given angle with respect to the normal to 
the sphere. This procedure is equivalent to injection of 
primary photons in different directions at fixed detector 
but is more timc-efhcient. For the sake of simplicity, we 
restrict our calculations here to the case when the source 
emits 7-rays isotropically. In this case the effects of dis- 
tortion of the circular symmetry of the surface brightness 
of the extended source due to the finite opening angle of 
the blazar jet and misalignment of the jet with respect to 
the line of sight do not affect the measurement of EGMF 
strength. We leave the investigation of the distortions of 
the morphology of extended emission caused by specific 
geometries of the jet for the future investigation. 

The results of numerical calculation of the energy- 
dependent morphology of extended emission arc shown 
in Figs. [4] and [5] for a source at a distance D = 120 Mpc 
(equal to the distance of the blazar Mrk 421). 

In the case of calculation shown in Fig. [4] the EGMF 
strength is assumed to be S = 10"^^ G. One can see that 
strong deflections of the cascade electrons by the EGMF 
lead in this case to a very large size of the extended emis- 
sion at the energies ~ 0.1 TeV, with significant contribu- 
tion to the signal extending up to ~ 4° from the source. 
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FIG. 4: The arrival directions of the primary and secondary 
cascade 7-rays (circles) from a source at a distance D — 
120 Mpc. The EGMF strength is 10"" G. The sizes of the 
circles representing each photon are proportional to the pho- 
ton energies. The blue dashed circle has radius 1.5°, equal to 
the radius of the FoV or MAGIC telescope. The radius of the 
blue solid circle is 2.5°, which corresponds to the size of the 
FoV of HESS telescope. 



Obviously, the extended emission at these energies can 
not be detected by a typical Cherenkov telescope array, 
which has a FoV with the diameter ~ 5° (HESS) or - 3° 
(MAGIC) (shown as, respectively, solid and dashed cir- 
cles in Fig. HI). The extension of the som-cc becomes 
smaller than the size of HESS and MAGIC FoVs only at 
the energies above ~ 1 TeV. 

The extended emission is more compact if the strength 
of EGMF is weaker. This is illustrated in Fig. [5] 
which shows the result of calculation of extended emis- 
sion around the same point source as in Fig. 31 but as- 
suming the EGMF strength B = 10~^^ G. One can see 
that in this case most of the extended emission down to 
the energies ~ 0.1 TeV fits into the FoVs of both HESS 
and MAGIC telescopes. 

Fig. [HI shows the radial surface brightness profiles of 
the extended sources shown in Figs. [4] and [5] in different 
energy bands. One can see that in each energy band the 
profiles of the extended emission produced by the cascade 
in a weaker magnetic field (open diamonds) are always 
more concentrated toward the central source than the 
profiles of extended emission produced by a cascade in 
a stronger magnetic field (filled diamonds, the magnetic 
field is 10~^^ G. The details of the shape of the surface 
brightness profile at the angular scale 9 ~ 0.1° ~ ^psf 
close to size of the PSF of a telescope could be calculated 
only if the details of the response of the telescope are 
known. 




FIG. 5: Same as Fig. Hbut for the EGMF strength 10""G. 



In our calculations we do not model the telescope's 
PSF and, instead, assume that 7-rays from a point source 
are homogeneously distributed within a disk of the radius 
^PSF = 0.1°. This explains the appearance of "plateau" 
in the surface brightness profiles. Comparing the nor- 
malization of the plateau with the one of the extended 
emission, one can judge, if the sensitivity of a telescope is 
sufficient for the detection of extended emission around 
a point source at a given level of the significance of de- 
tection of the source for a particular telescope. It is im- 
portant to note that, in spite of the much lower nor- 
malization of the surface brightness profiles of extended 
emission, as compared to the point source emission, the 
total flux of the extended source (integrated over the tele- 
scope's field of view) could be comparable to the one of 
the point source. 

The qualitative discussion of the model radial bright- 
ness profile of the extended emission, presented in the 
previous section, adopted simplifying assumptions (such 
as e.g. the monochromatic spectrum of the primary 7-ray 
source and monochromatic approximation for the spec- 
trum of IC emission from the cascade e~^e~ pairs). Al- 
though useful for the qualitative understanding of the 
method of measurement of EGMF, these assumptions do 
not hold in realistic situation. In a more general case, the 
shape of the radial surface brightness profile of extended 
source, S{6,E), could be derived from the results of our 
numerical calculations. 

In the particular case of a powerlaw primary 7-ray 
spectrum with the photon index F = 2, which we consider 
as an example, the radial brightness profiles of extended 
source shown in Fig. [5] can be approximated as powerlaw 
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FIG. 6: Surface brightness profiles of emission in different 
energy bands (I - O.fO-0.12 TeV, 2 - 0.12-0.f8 TeV, 3 - 
0.i8-0.27 TeV, 4 - 0.27-0.40 TeV) produced by tiie cascades 
in EGMF of tlie strengtfi fO^^G (thin lines) and fO^^^G 
EGMF (thick lines) and without magnetic field (red lines). 



with an exponential cutoff: 



dn^(E~f) 



oc 93(9, E-y) oc 



Bcut(£'7 



exp 



0cut(£'7 

(17) 

with the slope a ~ — 1. The energy-dependent param- 
eter 0cut(£'7) characterizes the extension of the source 
at a given energy E. This parameter is an analog of the 
parameter 0ext, which would characterize the size of the 
extended source in the case of the monochromatic pri- 
mary source spectrum. 



F. Measurement of EGMF strength 

Since the parameter Qcut, which appears in Eq. (|17p . 
is a direct analog of the parameter Qcxt for the case of a 
source emitting a powcrlaw spectrum of primary 7-rays, 
one expects that, similarly to Ooxt, ©cut should decrease 
inversely proportionally to the energy. Fig. [7] shows that 
this is indeed the case. 

The energy dependence of the size of the extended 
source, 9cut(^7) is determined by the strength of EGMF 
along the line of sight. Therefore, measuring this en- 
ergy dependence, one can extract information about the 
properties of EGMF from the observational data. In the 
simplest case of a monochromatic point source of 7-rays, 
the dependence of the source extension on the EGMF 
strength could be found analytically ©. It is clear that if 
the primary 7-rays from the source are distributed within 
a certain energy range, a superposition of the extended 
sources produced in result of absorption of the primary 
7-rays of different energies will lead to the changes of the 
normalization coefficient N = N{T,Ecut) depending on 
the spectral model of the primary source. Then Eq. ^ 



can be rewritten for the case of power law injection spec- 
trum via substitution 0oxt — ^ ©cut and introduction of 
Af(r,^cut): 



'dcutiEj, B) 



A^(r,^cut) 

t{E^„{E^),z) 



E. 



7 



6 TcV 



B 



10-13 G 



(18) 

Using a numerically calculated N{T, Ecut), one would be 
able to derive the strength of EGMF from the observa- 
tionally measured Qcut{E-y). 

To illustrate the implementation of the algorithm de- 
scribed above, we consider the case of a primary 7-ray 
source with the photon index F = 2 and and exponential 
high-energy cut-off at Ecut ~ 300 TeV. We derive the 
value of the normalization coefficient 7V(r = 2, Ecut = 
300 TeV) from a fit of the energy dependence of 8cut(-B-y) 
using the simulated properties of the extended emission 
for a (arbitrarily chosen) reference EGMF B = 10~^^ G. 
Once the value of the normalization coefficient N is fixed, 
one can resolve the Eq. (fT8|) with respect to B and use the 
measurements of ©cut (£'7) to infer the value of B from 
the simulated data sets produced for different choices of 
EGMF strength. 

The result of such a procedure of "reconstruction" of 
B from the simulated data sets is shown in Fig. [S] 
The y axis of the plot shows the reconstructed value of 
the EGMF strength, while the initially assumed EGMF 
strength is plotted along the x axis. Since the normal- 
ization coefficient N is found from the simulations as- 
suming the "reference" magnetic field, the reconstructed 
magnetic field is exactly equal to the assumed magnetic 
field at the reference point B = 1Q~^^ G. 

From Fig. [5] one can see that the assumed value of 
EGMF can be reconstructed from the data in a satis- 
factory way only in a certain range of EGMF strengths 
(10-1^ G < B < 10-12 G). This range is determined 
by the assumed characteristics of the 7-ray telescopes 
and/or by the intrinsic properties of the cascade. 

When the field strength gets close to 10"!^ G, the size 
of the extended source at the energies above ~ 0.1 TeV 
approaches the assumed angular resolution of the tele- 
scope, 9psp = 0.1°. Under these conditions, the extended 
emission becomes indistinguishable from the point source 
and measurement of ©cut(-E') becomes impossible. This 
is refiected in a large errorbar of the point at lO^^^ G in 
Fig. [5] Obvious ways to extend the range of accessible 
EGMF strengths toward still weaker fields is to improve 
the angular resolution and/or to extend the observed en- 
ergy band to the energies below 0.1 TeV. Another possi- 
bility is to observe such field with Fermi at the energies 
1-100 GeV. 

If the EGMF along the line of sight is stronger than 
lO^i'' G, the size of extended source at the energy 
^ 0.1 TeV becomes comparable to the size of the FoV 
of the current generation Cherenkov telescope arrays, 
^FoV ~ 2°. Under these conditions, the measurement 
of Qcut{E) is still partially possible at the energies much 
larger than ^0.1 TeV, but the precision of determination 
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of Qcut{E) rapidly drops when the field becomes larger 
than ~ 10~^^ G. One should note that, in a realistic situ- 
ation, the measurement of the radial brightness profile of 
extended source becomes problematic already when the 
source size is smaller, but comparable to the size of the 
field of view. The difficulty is that the source signal has 
to be disentangled from the cosmic ray background, in 
the camera, with intensity depends on the off-axis angle. 

In our calculations we have ignored these difficulties 
and assumed that (a) the telescope FoV is not limited 
and (b) that an arbitrary large source (with extension 
up to ~ 10°) could be detected. One can see that with 
a large FoV, high sensitivity telescope the field strength 
can be reconstructed from the measurement of character- 
istics of the extended emission up to the EGMF strengths 
^ 10"^2 G, at which the energy dependence of the source 
size on EGMF gets lost because of the too strong deflec- 
tions of the cascade electrons and positrons. 

In principle, an energy-dependent size of the extended 
source is predicted also when the strength of EGMF is 
larger than 10~^^ G, within the model of extended emis- 
sion discussed in the Ref. 0. However, in this case 
the size of the source is determined by the mean free 
path of the 10-100 TeV 7-rays through the cosmologi- 
cal infrared/optical background. The angular extension 
of the source in this model can be readily estimated to 
be Q{Ej) ~ [T{E^g)]~^. Taking into account the fact 
that in a limited energy range T{E^g) could be approxi- 
mated by a linear function, one finds that the size of the 
extended source is also expected to shrink with energy. 
Contrary to the case of the EGMF-dependent source size, 
the size of the halos considered in the Ref. decreases 
as e{E^) ~ E-^ ~ E~^^^ (see Eq. @). This means 
that the energy-dependence of the size of extended halos 
around extragalactic sources considered in the Ref. 01 
could be readily distinguished from the one expected if 
the extended source size is determined by the EGMF. 
Systematic detection of the extended emission with the 

— 1/2 

characteristic E-y dependence of the source extension 
on the energy would provide a strong argument in favor 
of the existence of strong EGMFs, B > 10"^^ G. 

At the end of this section, we would like to comment 
on the implementation of the proposed method for the 
analysis of the data of the real 7-ray observations of ex- 
tragalactic sources. In a realistic situation the primary 
source spectrum is (a) not known in advance and (b) 
is strongly modified at the highest energies by the ef- 
fect of absorption on the EBL. This means that, in gen- 
eral, one can not assume a particular spectral model of 
the primary sources for the EGMF search analysis, as it 
is done above for the analysis of the Monte-Carlo sim- 
ulated data sets. Instead, the measurement of EGMF 
has to be done simultaneously with the measurement of 
the parameters of the high-energy emission spectrum of 
the primary source. This implies the following analy- 
sis procedure. First, using the Montc-Carlo simulations, 
described above, one has to produce a "look-up table" 
of the normalization coefficients N{T, Ecut, ■■■) for a set 




0.1 1 
E^TeV 

FIG. 7: Gcut as function of energy E.y for several values of 
EGMF, lO"^'^ G (squares), 2-10-^^ (circles), 5- 10"^^ (trian- 
gles), 10~^* (diamonds). Solid lines correspond to the fits of 
the data points with Ocut ~ E^^ . 




FIG. 8: Comparison of the EGMF strength reconstructed 
from measurement of extended emission in simulated data 
sets, -Brcconstr, with the real value of the magnetic field 
strength assumed for the simulations, -Breai- The line cor- 
responds to Broal = -Brcconstr 



of parameters (F, i?cut, •■•■) of different possible spectral 
models of the 7-ray source. Next, using these look-up 
tables, one can simultaneously find the best-fit values 
for the spectral model parameters (T, Ecut, ■■■) and for 
the EGMF strength B, via a multi-parameter fitting 
the spectral and morphological characteristics of the ex- 
tended emission. 
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IV. SUMMARY AND CONCLUSIONS 

111 this paper we have used Monte-Carlo simulations 
of development of 7-ray induced electromagnetic cascade 
in the intergalactic space to study the properties of ex- 
tended emission around extragalactic point sources of 
VHE 7-rays. We have demonstrated that the proper- 
tics of the cascade and, as a consequence, the properties 
of the expected extended emission around extragalac- 
tic 7-ray sources strongly depend on the strength of the 
(largely uncertain) EGMF along the line of sight toward 
the source. 

The dependence of the properties of the cascades on 
the EGMF strength could be used to reveal the pres- 
ence of extremely weak magnetic fields in the intergalac- 
tic medium and to measure (or, at least, estimate) the 
strength of these fields. Such a measurement would en- 
able to distinguish between existing competing cosmo- 
logical models of the origin of magnetic fields in the Uni- 
verse, which widely differ (by more than 10 orders of 
magnitude) in their predictions of the EGMF strength. 

We have demonstrated that if the typical magnetic 
field strength along the line of sight is in the range 
10~^^ G< B < 10~^^ G, the cascade emission could be 
detected as an energy-dependent extended " glow" around 
the point source (see Figs. SI O by the ground-based 7- 
ray telescopes. Telescopes optimized for the searches of 
stronger and weaker magnetic fields should have quite dif- 
ferent characteristics. Large effective collection area and 
large size of the telescope's field of view are needed for the 
detection of stronger magnetic fields B 10"^^ G. At the 
same time, high sensitivity at low 7-ray energies (below 
0.1 TeV) and good angular resolution are important for 
the detection of weaker B ~ 10^^^ G fields. Still weaker 
fields B < 10~^^ G should reveal themselves through the 
extended emission around extragalactic 7-ray sources in 
the Fermi energy band, ^ 1 — 100 GcV. 

The extension of the 7-ray glow around extragalactic 
point sources appears to decrease inversely proportion- 
ally to the energy (Figs. [SI E]). Measurement of the 
energy dependence of the source size enables to derive a 



constraint on the strength of extragalactic magnetic field 
from the data (Monte-Carlo simulated data in our case), 
using Eq. (fTT]) . (fT8|) . We have described the algorithm of 
reconstruction of the EGMF strength from the data. Fig. 
|8] presents the results of implementation of the algorithm 
for the analysis of the simulated data sets. The proposed 
algorithm could be directly used for the analysis of the 
real data of ground and space- based 7-ray telescopes. 

Apart form the measurement of the EGMF, a study of 
the properties of extended emission around extragalac- 
tic point sources enables to constrain the properties of 
the primary 7-ray source spectrum at the energies above 
10 TeV at which the flux from the source is strongly sup- 
pressed by absorption on the extragalactic background 
light. We have argued that, in fact, in the analysis of the 
real data, the constraints on the spectral characteristics 
of the primary source have to be found simultaneously 
with the constraints on the EGMF strength. 

When this paper was finished, we became aware of a 
paper which analyzes the properties of extended emission 
around extragalactic 7-ray sources by Dolag et al. [36j . 
which presents results complementary to our study. Con- 
trary to our study, the Ref. [3^ concentrates on the anal- 
ysis of the properties of the EGMF-dependent extended 
emission assuming different spectral and morphological 
characteristics of the primary 7-ray source, fixing a par- 
ticular model of EGMF. 
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